Vortex matter in the charged Bose liquid at absolute zero 



o 
o 

O 

Q 



V. V. Kabanov 1 and A. S. Alexandrov 2 
1 Josef Stefan Institute 1001, Ljubljana, Slovenia 
2 Department of Physics, Loughborough University, Loughborough, United Kingdom 

The Gross-Pitaevskii-type equation is solved for the charge Bose liquid in an external magnetic 
field at zero temperature. There is a vortex lattice with locally broken charge neutrality. Remark- 
ably, there is no upper critical field at zero temperature, so the density of single flux-quantum 
vortices monotonously increases with the magnetic field up to B — oo and no indication of a phase 
transition. The size of each vortex core decreases as about B -1 ' 2 keeping the system globally charge 
neutral. If bosons are composed of two fermions, a phase transition to a spin-polarized Fermi liquid 
at some magnetic field larger than the pair-breaking field is predicted. 

PACS numbers: 74.20.-z, 74.65. +n,74.60.Mj 
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Charged Bose liquids (CBLs) have been solely of aca- 
demic interest for a long time 

Notwithstanding, experimental realization of the Bosc- 
Einstein condensation (BEC) of trapped ultra-cold atoms 
El El H3. H3 . Ilif made it possible to create ultra- 
cold plasmas |15| by using lasers to trap and cool neu- 
tral atoms to temperatures of 1 mK or lower. Another 
laser then ionizes the atoms by giving each of the out- 
ermost electrons just enough energy to escape the elec- 
trical attraction of its parent ion. The ions retain the 
millikelvin temperatures of the neutral atoms, so they 
may bose-condense, if their spin is an integer. There 
is also growing experimental evidence for preformed 2e- 
charged bosons in high-temperature cuprate supercon- 
ductors, such as normal state pseudogaps, unusual up- 
per critical fields, small normal-state Lorentz numbers, 
etc 0| . Similar charged boson physics is expected in 
a lattice of mesoscale superconducting dots, if parame- 
ters are chosen such that single-electron tunnelling is sup- 
pressed and only Cooper-pair char ges tunnel between the 
domains via Josephson tunnelling • It is also possible 
to describe the universal features of the superconductor- 
insulator transition as a function of disorder in quasi-two 
dimensional systems in terms of boson physics [Ta . ll9j . In 
order to model the transition in terms of bosons, one has 
to include the Coulomb repulsion, otherwise all bosons 
would collapse into the lowest lying highly localised state. 

These developments have renewed interest in CBL as a 
fundamental reference system. A non-interacting gas of 
charged bosons cannot bose-condense at any finite mag- 
netic field because of a one-dimensional character of mo- 
tion in the lowest Landau band pj. However, interact- 
ing charged bosons condense below some (upper) critical 
field B ^ H C 2(T) since their collisions remove the one- 
dimensional singularity of the density of states l20l. The 
BEC field diverges with decreasing temperature |20L l2ll| , 
so that H C 2{T) — oo at absolute zero. A single vortex 
in CBL has a charged core and an electric field inside 
[23 , while its ma gnet ic field is virtually identical to the 
Abrikosov vortex [23 . 

Here we present the ground state of CBL in an ar- 
bitrary magnetic field solving numerically the Gross- 
Pitaevskii -type equations with the long-range Coulomb 



interaction between bosons. We find a lattice of charged 
vortices, which does not disappear at any finite magnetic 
field. The density of vortices monotonously increases 
and their core size decreases with the magnetic field up 
to B = 00. The size of vortices also depends on the 
thickness of CBL films different from the conventional su- 
perconducting films. When bosons are composed of two 
fermions, there is a phase transition to a spin-polarized 
Fermi liquid at some magnetic field. 

The Hamiltonian of charged bosons on a compensat- 
ing homogeneous background (to ensure global charge 
neutrality) in the external magnetic field with the vector 
potential A(r) is given by 
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where m,e,n,fi are the mass, charge, average density and 
chemical potential of bosons, respectively, and V(r) = 
e 2 /|r| is their Coulomb repulsion |24| . 

The equation of motion for the Heisenberg field opera- 
tor, ip(r, t), is derived using this Hamiltonian. If the den- 
sity is relatively high, so that the dimensionless Coulomb 
repulsion r s = me 2 //i 2 (47rn/3) 1 / 3 is not large, one can ex- 
pect that the occupation numbers of one-particle states 
are not very much different from those in the ideal Bosc- 
gas. In particular, one state remains macroscopi cally oc- 
cupied at T — OK. Then, following Bogoliubov |23 one 
separates the large matrix element from tp by treating 
the rest yb as small fluctuations, i/j(r, t) = ipo(r) + ip(r, t). 
The anomalous average ?/>o( r ) = (^( r >*)) is approxi- 
mately equal to ^/n in a homogeneous system at T = 0. 
Substituting the Bogoliubov displacement transforma- 
tion into the equation of motion and collecting c— number 
terms of i/jq, one obtains the equation for the macroscopic 
condensate wave function as [22] 
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FIG. 1: A few vortices in a sample of the size L = 22. 4^, (a). 
The phase profile of the order parameter is shown in (b) . The 
phase changes by 2n around each core at any magnetic field. 




FIG. 2: Single vortex in CBL [22| (a) compared with the 
Abrikosov vortex (b). 
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The integro-differential equation (2) is quite different 
from original Ginzburg-Landau fGLl[26| and the Gross- 
Pitaevskii |27| equations, describing the order parame- 
ter in the BCS and neutral superfluids, respectively. As 
recognised by one of us the coherence length in CBL 
is just the same as the screening radius, so the core 
of a single vortex is charged [2^. Indeed, introducing 
the dimensionless quantities: / = ^^/n 1 / 2 , p = r/A, 
h = e^XcurlA/hc for the order parameter, length and 
magnetic field, respectively, one obtains the coherence 
length about the same as the screening radius at T = OK 
[I, £ = (fi/2 1 / 2 m*w p ) 1 / 2 , where lo p = Unne 2 /m) 1 ' 2 is 
the zero-temperature plasma frequency [2j. The London 
penetration depth is conventional, A = (mc 2 /47rne 2 ) 1 ^ 2 , 
but a new feature is an electric field potential, e<p(r) = 
J dr'V(r — r')[|-0o( r ')| 2 — n]. Moreover, the chemical po- 
tential fx is zero, as it should be in the globally neutral 
CBL in the thermal equilibrium below the BEC critical 
temperature. 

Any realistic CBL is an extreme type II with a very 
large Ginzburg-Landau parameter, n = A/£ 1 [22j . 
For example, the coherence length and the electric field 
inside the vortex core are about 1 nm or less and 10 mV, 
respectively, with the material parameters typical for 
cuprates (m = 10m e , n — 10 21 cm~ 3 and eo > 100), and n 
is about 10 3 with these parameters. Hence, the magnetic 
field is practically homogeneous, and the ground state 
■00 ( r ) can be found by minimizing the energy functional 
E with respect to ipo(r), 

dr [ dr'^(r-r')|^o(r)| 2 (|Vo(r')| 2 -2n),(3) 



where A = {0, Bx,0}. In numerical simulations we con- 
sider a sample with the rectangular cross-section L x L 
and the magnetic flux BL 2 = m<&o, where m is an inte- 
ger ($ is the flux quantum). When the magnetic field B 
is applied along z-direction, the order parameter if>o{x, y) 
does not depend on z obeying the following translation 



symmetry, 

ipo(x +L,y) = exp(-ieBLy/tlc)ipo(x,y) 
tpo(x,y + L) = ip (x,y). (4) 

These relations can be used as boundary conditions when 
m is an integer. 

Because Eq.(3) does not contain the penetration 
depth, it is convenient to introduce new dimension- 
less coordinates x = r/£, the vector potential, a = 
(0, 27rf?2;£ 2 /$o, 0), and the Coulomb energy u(x) = 
e</>/(w 2 m*£ 2 ). As a result, the problem is reduced to 
minimization of the functional 
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where the Coulomb field satisfies the Poisson equation, 

At-(x) = l-|/(x)| 2 . (6) 

To compare CBL vortex state with the Abrikosov vortex 
lattice we also minimize the conventional GL functional 
using the same dimensionless unites, 
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where £ = h 2 /(2m\a\) 1 ^ 2 , n s — \a\/(3 and the order pa- 
rameter / is normalised by^/nj. Here a and (3 are con- 
ventional GL coefficients |2g. We apply the standard 
discretization procedure described in Ref. 29]. Eq.(6) for 
the electrostatic potential is solved by the Fourier trans- 
form in the discrete form, and the resulting energy is 
minimized with the conjugated gradient algorithm. 

Since both functionals depend only on the dimension- 
less vector-potential a which is proportional to the prod- 
uct B^ 2 , simulations can be performed at fixed L and £ 
by changing B or at fixed L and B by changing £. Our 
numerical results are shown in Figs. 1-4. At any value of 
the magnetic field we find the triangular vortex lattice. 
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FIG. 3: The vortex lattice in CBL for 30 flux quantum 
per cross-section, (a) L/£ = 33.67, (b) L/£ = 25.25, (c) 
L/£ = 14.43, (d) L/£ = 10.1 One can see from the scale near 
each figure that the order parameter remains large outside 
the cores, / > 1 at any £ (or at any magnetic field). 



FIG. 4: The Abrikosov vortex lattice for 30 flux quantum per 
cross-section (a) L/£ = 33.67, (b) L/£ = 25.25, (c) L/£ = 
14.43, and for (c) L/£ = 13.87, which corresponds to B close 
to H C 2 = $o/(27r^ 2 ). The order parameter decreases when B 
approaches the conventional upper critical field. 



While the field is small, there are only a few vortices 
per sample cross-section, Fig. la. When vortices are far 
apart, their interaction yields only a small contribution 
to the energy functional but even in that case a triangular 
lattice of vortices is clearly seen in CBL, Fig.l. 

Each vortex carries one flux quantum, as can be seen 
from the phase profile in Fig. lb. It has an unusual core, 
Fig. 2, in agreement with Ref. j^, which differs qual- 
itatively from the Abrikosov vortex |2^| due to a local 
charge redistribution caused by the magnetic field. The 
breakdown of the local charge neutrality, Fig. 2, is due to 
the absence of an equilibrium normal state solution in 
CBL at T = with ^ = 0, as explained in Ref.^. 

Increasing the field first increases the vortex density 
with about constant size of the cores, as in conventional 
superconductors, Fig. 3 and Fig. 4. However, quite differ- 
ent from the Abrikosov lattice, increasing the field further 
does not lead to a superfluid to normal phase transition, 
but instead it increases the density of vortices by de- 
creasing the size of every individual core, Fig. 3c, d. Re- 
markably, each vortex carries one flux quantum at any 
field. Keeping the global charge neutrality the charge 
heterogeneity depends on the magnetic field, and the 
core diameters decrease with the field, when the field 
is large, £ 2 > 2ttHc/ (eB). Indeed, in this regime the 
"bare" coherence length £ becomes irrelevant, but the 
only characteristic length is the distance between single 
flux-quantum vortices, i.e. r w yj2nhc/ {eB). As a re- 
sult, the amplitude real-space modulations of the order 
parameter increase with the magnetic field in CBL, while 
they decrease in conventional superconductors, where 
the order parameter vanishes at and above the finite 
H c2 = $ /(2< 2 ) (fig.4c,d). 

There is another difference between CBL and conven- 



tional vortex matter in case of thin films. If we assume 
that the film thickness d is small, d <§C £, then the left 
hand side of Eq.(6) takes the form (1 — \f(x.)\ 2 )dS(z)/^. 
The dimension analysis readily shows that the true co- 
herence length, ^2D depends on the thickness as = 
(^/(i) 1 / 3 in that case. As a result the size of vortex cores 
depends on the thickness of CBL films different from the 
conventional films. 

There is also an important consequence of the infi- 
nite (orbital) upper critical field at absolute zero in such 
CBLs, where singlet bosons are formed of two fermions 
In this case sufficiently large magnetic field can 
break bound pairs via a spin-flip of one of two fermions, 
if triplets are unstable, because the singlet binding en- 
ergy A decreases with the field as A(B) = A — 2ugB 
(jiB = eh/(2m e ) is the electron Bohr magneton) |30| . 
A spin-polarised Fermi liquid appears at B ^ H p , where 
Hp — A/(2/ie) is the pair-breaking field. In this estimate 
we neglect the orbital (Landau) diamagnetism of bosons 
and fermions, and the Coulomb energy of the charged- 
modulated vortex lattice. The latter is of the order of 
ecf> c nt; 2 B/$o per unit volume, where <fi c ~ fi, 2 /(em£ 2 ) 
is the characteristic electrostatic potential inside vortex 
cores. The Coulomb energy is small compared with the 
spin (Pauli) contribution if m e /m <C 1, which we as- 
sume to be the case, so diamagnetic contributions are 
also small. However, bound pairs still survive up to a 
higher field H* = H p + ti/(Nhb) > H v due to the Pauli 
exclusion principle, which prevents any further decay of 
pairs, if the number of fermions « N\A(B)\ remains 
smaller than 2n ( N is the fermion density of states). 
There is a boson- fermion mixture, if H p < B < H*, with 
the fermion density modulated in real space because of 
charged vortices. Normal fermions (as well as normal 
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bosons pushed up from the condensate by temperature) 
are distributed inhomogeneously across the sample with 
the maximum density in the vortex cores, where their 
potential energy is at minimum. The excess density of 
normal carriers inside the cores screens the electric field 
caused by the inhomogeneous condensate density. If the 
screening length due to normal fermions becomes smaller 
than the coherence length £, one can expect a nontrivial 
field dependence of the size of vortices, which disappear 
at B=H*. 

In conclusion, we have found the triangular lattice of 
single-flux-quantum charged vortices in CBL which can- 
not be destroyed by any magnetic field at zero temper- 
ature. The vortex density monotonously increases and 
their core size decreases with the magnetic field up to 
B = oo with no indication of a phase transition. The 
core size depends on the thickness of CBL films. At 
finite temperatures H C 2 (T) is finite [2(], H]] . Neverthe- 



less, unusually large charge modulations with the scale 
depending on the magnetic field should persist at finite 
temperatures as well. The phase transition to the spin- 
polarized Fermi liquid at some magnetic field larger than 
the pair-breaking field has been predicted for preformed 
bosonic pairs. These results are relevant for real charged 
Bose-liguids in ultracold plasmas and in the supercon- 
ducting cuprates, and for superconducting quantum dots 
and superconductor-insulator phase transitions described 
by a similar boson physics. There is also a close analogy 
between the vortex structure in CBL and the Josephson 
vortices. Since the normal phase is not defined below T c , 
there is no "normal" vortex core in CBL, and there is 
no "normal" core in the Josephson vortex either. One 
can define the lower critical field H c \ when a first vortex 
penetrates into CBL |22| and into the Josephson junction 
[3l|, but the upper critical field is infinite in both cases. 
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